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We report test results of the correlation between time variability and peak luminosity of 
Gamma-Ray Bursts (GRBs), using a larger sample (32) of GRBs with known redshift than that 
available to Reichart et al. (2001), and using as variability measure that introduced by these 
authors. The results are puzzling. Assuming an isotropic-equivalent peak luminosity, as done 
by Reichart et al. (2001), a correlation is still found, but it is less relevant, and inconsistent 
with a power law as previously reported. Assuming as peak luminosity that corrected for GRB 
beaming for a subset of 16 GRBs with known beaming angle, the correlation becomes little 
less significant. 
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1 INTRODUCTION 

Despite the small number of Gamma-Ray Bursts (GRBs) with 
known redshift (several dozens), several correlations between in- 
trinsic temporal or spectral parameters of the GRB prompt emission 
and GRB energetics have been discovered in the last seven years. 
[NoitIs et al. ( 2000) found an anticorrelation between peak luminos- 
ity and the spectral lag (obtained by cross-correlating the time pro- 
files of the same GRB in various energy ban ds), according to which 
more luminous bursts exhibit shorter lags. ISalmonson & GalamJ 
i2002 h discovered a positive correlation between the spectral lag of 
the gamma-ray prompt emission and the jet-break time of the after- 
glow decay, according to which a small break time corresponds to 
a small lag and consequently to a high peak luminosity of the GRB. 
Concerning the temporal properties of GRB time profiles, evidence 
has been found for a positive correlation between temporal vari- 
ability of the light curves and isotropic-equivalent peak luminosity 
for the GRBs with known redshift (Reichart et al. 2001, hereafter 
ROl; Fenimore & Ramirez-Ruiz 2000). 

Moreover Reichart et al. (2003) have shown that the variabil- 
ity vs. peak luminosity correlation could also hold true for X-Ray 
Flashes (XRFs; see Heise et al. (2001)). As a consequence of the 
mentioned correlations, a correlation between time variability and 
spectral lag is also expected and confirmed for a large sample of 
BATSE bursts ( Schaefer et al. 2001). The variability vs. peak lu- 
minosity correlation has been explained by several authors (e.g., 
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iKobavashi etaP te002'):'M eszaros et all <2002h ) mainly within the 
framework of the standard fireball model, according to which inter- 
nal shocks between ultra relativistic shells are responsible for the 
pulse-like structure of the GRB prompt emission, while the smooth 
afterglow emission is due to external shocks between the fi reball 
wind and the matter surrounding the GRB progenitor fe.g.. |PirM] 
>2004) for a review). 

GRB variability-connected properties are thought to be more 
sensitive to the bulk Lorentz factor V and, if the GRB emission is 
beamed, to the jet ope ning a n gle and/or the viewing an gle (e.g., 
ISalmonson & Galama; <20(32h : iloka & Nakamur j <200ll) ). Within 
the fireball model, there are different mechanisms that could ac- 
count for different time variabilities, also giving possible explana- 
tions for XRF properties ( Meszaros et al. 2002). In addition, the 
"cannonball model" for GRBs ( «Dado. Dar & De Ruiula.20 02) also 
seems to ex plain the variability vs. peak luminosity correlation 
<Plagj200ll) . 

From the above correlations als o luminosity e stima- 
tors have tentatively been derived jReich art et alj 1200 it 
'Fenimore & Ramirez-Rui3 l200(t ISchaefer et alj 1200 ll) to in- 
vestigate general properties of GRBs, such as the luminosity 
function and the possible link with the star formation rate. In 
addition, empirical redshift indicators have been proposed based 
on the calibration derived with the small sample of GRBs with 
known redshift, making use of the X- and gamma-ray observations 
alone ( Atteia 2003; Bagolv et al. 2003). 

In this work we test the variability vs. peak luminosity cor- 
relation using the variability definition given by ROl. We used a 
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sample of 32 GRBs with known redshift. Furthermore we studied 
the same correlation by replacing the isotropic-equivalent peak lu- 
minosity with that corrected for beaming for a subset of 16 GRBs 
with known collimation angle provided by Ghirlanda et al. ( 2004). 

In Section|2| we discuss our sample of GRBs; in Section|3]we 
discuss the time variability analysis; in Section |4| we estimate the 
peak luminosity of the GRB in our sample and compare it with the 
ROl results. In Section|5|we present our results on variability/peak 
luminosity correlation, in Section|6|we discuss our results. 



2 THE GRB SAMPLE 

2.1 GRBs with known redshift 

The sample of 32 GRBs with known redshift includ es 16 GRBs 
detec t ed by the Gamma-Ray Burst Monitor (G RBM) jperoci et alJ 
1 1 9971: |Frpnteraetal."l997"; 'Costa et al."l998') on board the Bep- 
poSAX satellite iBoellaetal. 199 7) during the period 1997-2002, 
two by the BATSE experiment jPaciesas et alJll999t) aboard the 
Compton Gamma-Ray Observatory (CGRO), six by the FRE- 
GATE instru ment aboard HETE- II jAtteia et alj l2003h. one by 
Konus/WIND ( ADtekar et al. 1995), one by Ulysses iHurlev et al] 
[l992) and six by BAT/ Swift i Gehrels 2004J . Eight of the 16 
GRBs detected with the BeppoSAX GRBM were also detected 
with BATSE. We used public archives for GRB data obtained with 
BATSE\ HETE-lf, Konus/WIND'-^, and BAT/Swi/f^. TableQre- 
ports the list of the GRBs in our sample with mentioned the space- 
craft that detected it. When the same GRB has been detected by 
more than one instrument, we first checked the consistency of the 
results derived from different data sets and then concentrated on the 
instrument which had the best SNR. 

The time binning of the GRB light curves in our sample, 
which was used to derive the time variability, was the following: 
7.8125 ms for the GRBM data, 64 ms for BATSE, 164 ms for 
HETE-II, 64 ms for Konus/WIND, 31.25 ms for Ulysses. In the 
case of BAT/Swift we made use of the event files and extracted 
the mask-tagged light curves with a binning time of 8 ms. Given 
that for the GRBs detected with the BeppoSAX GRBM, the high- 
resolution (7.8125-ms binning) time profiles are available in the 
40-700 keV energy band, for the others we used the light curves in 
the energy bands which have the largest intersection with the 40- 
700 keV band: 1 10-320 keV (channel 3) for BATSE, 30-400 keV 
(band C) for FREGATE, 25-100 keV for Ulysses and 50-200 keV 
for Konus/WIND. In order to match the GRBM band, for BAT/Swift 
we extracted the light curves from the event files in the 40-350 keV 
band. 

For GRB990510, given that the 7.8125-ms GRBM light curve 
is not available, we preferred to use 64-ms BATSE data rather than 
the 1-s GRBM light curve. 

Six GRBs (980613, 011211, 021004, 050126, 050318 and 
050416) with known redshift were not included in our sample due 
to their low signal, which prevented us from deriving a statistically 
significant variability estimate. Another GRB (021211, detected 
with HETE-II) was not included in the sample, due to the high ra- 
tio between binning time and smoothing time which could bias the 

^ ftp://cossc.gsfc.nasa.gov/compton/data/batse/ascii_data/64ms/ 

2 http://space.mit.edu/HETE/Bursts/Data/ 

^ http://lheawww.gsfc.nasa.gov/docs/gamcosray/legr/ 

bacodine/konus.grbs.html 

^ http://swift.gsfc.nasa.gov/ 



variability estimate (see Section|3j. Other GRBs with known red- 
shift (00030IC, 000418, 000926) detected by both Konus/WIND 
and Ulysses were not included in our sample, because unfortunately 
both Konus public and Ulysses data cover their light curves only 
partially. 

In the case of BATSE (970828 and 000131) the usual 4- 
channel 64-ms light curves are not available. Thus we made use 
of the 16-channel MER spectra acquired along either entire GRB 
with an integration time of 64 ms. Therefore we rebinned the 16- 
energy-channel MER data both spectrally and temporally in order 
to reproduce as much as possible the 4-channel scheme of BATSE 
64-ms light profiles. As we discuss below, we relied on coarse time 
resolution light curves only when the overall duration of the GRB 
was very long compared to the binning time. 

In general the data available cover the entire time profile of 
the GRBs in our sample. However there are some exceptions. In 
the case of the GRBM events, given that the high-resolution data 
cover 8 s before the trigger time and 98 s after it, in the case of 
the longest events (990506 and 010222, Tgo = 129 s and Tgo = 
97 s, respectively), it is not true. In these cases, the measure of 
time variability was obtained summing the variability in the part 
covered by the 7.8125-ms bins with that in the part covered by 1-s 
ratemeters (the tail of the burst). 

GRB000210 (Tgo = 8.1 s) suffers from a 2.5-s long gap due 
to corrupted high-resolution data occurred in the middle of the burst 
profile. Using the 1-s data, the mean 7.8125-ms light profile in the 
gap was reconstructed adding to the mean profile a Poisson noise. 
The value of GRB time variability so derived was found not to sig- 
nificantly change, even adding in the gap a non-Poisson noise com- 
patible with the 1-s time profile. 

For each GRB detected with the GRBM, we considered the 
light curves of the two most illuminated units and checking whether 
the best signal-to-noise ratio (SNR) was obtained from a single unit 
or by summing the two units. 

We found that for t he 11 GR Bs detected by GRBM and Wide 
Field Cameras (WFCs) ijager et al..l997.) . the best signal is ob- 
tained from a single GRBM unit (that with the larger area exposed 
to the GRB). For the five bursts detected with the GRBM but not 
with the WFCs the best signal was obtained by summing the two 
most illuminated units: units 1 and 4 (980703), 3 and 4 (990506, 
020405), 2 and 4 (991216, 010921). In principle the operation of 
adding the counts of different units is questionable because of dead 
time, as it will be discussed in Section|3| In practice, for the above 
cases we made sure that the results were consistent with those ob- 
tained when considering only the most illuminated unit for each 
GRB. This has been found to be no longer true, i.e. the correction 
for dead time becomes not negligible, when considering very small 
smoothing time-scales (Rossi et al., in preparation). 

As far as the 8 bursts detected with both GRBM and BATSE 
(970508, 971214, 980425, 980703, 990123, 990506, 990510, 
991216) are concerned, we used the BeppoSAX data for 971214, 
980703, 990123, 990506, and 991216, for which the higher time 
resolution of the GRBM turned out to be essential for a better vari- 
ability estimate, while for the remaining 3 GRBs (970508, 980425, 
and 990510) we used the BATSE data given the better SNR, after 
verifying the mutual consistency of the GRBM and BATSE vari- 
ability results. 
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Table 1. Variability vs. Peak Luminosity for 32 GRBs witli known redsliift (Icr errors). 



GRB 
Name 



Mission(^) 



Redshift 



T>=0.45 



Peak Lum. L'*') 
(10™ ergs-i) 



:Refs.('=) 



970228 


0.695 


BS/U/K 


2.2 





fj„„_|_0.018 
^^■^-0.017 


48.7 ± 9.9 


1 


970508 


0.835 


B/BS/U/K 


2.4 





„„„_Ln 013 


9.43 ± 1.89 


2 


970828 


0.958 


BAJ/K/S 


12.9 





-LQ 002 


120.0 ± 40.0 


3 


971214 


3.418 


BS/B/U/K/N/R 


4.4 





^ 012 
^^^^-0.012 


360. ± 65. 


4 


980425 


0.0085 


B/BS/U/K 


4.7 





„ 048 
^^^-0.048 


0.0007 ± 0.0002 


5 


980703 


0.966 


BS/B/U/K/R 


3.2 





„ . . -LQ 007 
^^^^-0.007 


26.4 ± 5.6 


6 


990123 


1.6 


BS/B/U/K 


12.8 





fj+0.002 


840. ± 121. 


7 


990506 


1.3 


BS/B/U/K/R 


8.6 





_ _^0.005 
^'^-0.005 


583. ± 121. 


8 


990510 


1.619 


B/BS/U/K/N 


3.2 





A -1-0 005 

^-L^-0.008 


300. ± 50. 


9 


990705 


0.86 


BS/U/K/N 


8.0 





+0.003 
^'°-0.003 


134. ± 21. 


10,11 


990712 


0.434 


BS/U/K 


4.1 





„ .r,-\-0 017 

^^^^-0.017 


5.4 ± 1.0 


12 


991208 


0.706 


K/U/N 


5.1 





_^0.003 
Uo^_0.003 


290. ± 100. 


13 


991216 


1.02 


BS/B/U/N 


2.6 





^o+O 002 
^^■^-0.002 


1398. ± 200. 


14 


000131 


4.5 


B/U/K/N 


8.0 





o„+0.005 

-■-^'-o.oos 


3600. ± 900. 


15 


000210 


0.846 


BS/U/K 


1.59 





„„„_Ln 002 
026+^™^ 


480. ± 50. 


16 


000911 


1.058 


U/K/N 


5.2 





^„+0.034 


360. ± 60. 


17 


010222 


1.477 


BS/U/K 


6.62 





_i_0 003 
201Io:oo3 


801. ± 119. 


18 


010921 


0.45 


BS/H/U/K 


5.3 





^„^_Ln 016 
"•^"-0.016 


8.0 ± 2.0 


19 


011121 


0.36 


BS/U/K/O 


8.3 





049+^S°^ 


19.9 ± 3.1 


20 


020124 


3.198 


H/U/K 


8.8 





,^•3+0.031 
^"■^-0.032 


300. ± 60. 


21 


020405 


0.69 


BS/U/K/O 


9.9 





ico+O.OOT 
^"°-0.007 


71.4 ± 11.2 


22 


020813 


1.25 


H/U/K/O 


17.4 





r,-o+0.007 
^^"-0.007 


340. ± 70. 


23 


030226 


1.98 


H/K/0 


26.6 





Q42+0015 
^^^-0.015 


25.0 ±5.0 


24 


030328 


1.52 


H/U/K 


24.9 





ncri +0.005 
"''^-O.OOS 


90. ± 18. 


25 


030329 


0.168 


H/U/K/O/RH 


4.9 





,(,,.+0.007 
^"^-0.007 


6.1 ± 1.2 


26 


041006 


0.712 


H/K/RH 


8.0 





nc9+0.002 
'^"^^-0.002 


66. ± 10. 


27 


050315 


1.949 


BSw 


12.3 





"^^-0.031 


38. ± 8. 


28 


050319 


3.24 


BSw 


3.6 





Qg-,+0.032 
-0.030 


84. ± 20. 


29 


050401 


2.90 


BSw 


4.4 





iqc;+0.028 
^^^^-0.029 


740. ± 100. 


30 


050505 


4.27 


BSw 


9.0 





■ppi ,.+0.043 
^^^-0.044 


250. ± 50. 


31 


050525 


0.606 


BSw 


2.0 





„ 1+0.003 
^^^-0.003 


80. ± 10. 


32 


050603 


2.821 


BSw 


1.2 





94^+0.037 
^^■^-0.034 


1200. ± 200. 


33 



t'') Mission: BS (BeppoSAX), B (BATSE/CGRO), K (Konus/WWD), H (HETE-II), U (Ulysses), S (SROSS-C), N (NEAR), R 
(RossiXTE), O (Mars Odyssey), RH (RHESSI), BSw (JiA^T /Swift): the data used are taken from the first mission mentioned. 
C') Isotropic-equivalent peak luminosity in 10^*^ erg s"'^ in the rest-frame 100-1000 keV band, for peak fluxes measured on a 1-s 
time-scale, Hq = 65 km s~^ Mpc~^, Sim = 0.3, and = 0-7. 

^'^^ References f or the redshift measurements: (1) Diorsovski et al. ( 1999"), (2)"Metzser et al. 1 1997), (3) Diorgovski et al. 1 2001^, 
r4)lKulkami et aL 1 1998) . (5) Tinnev et al. |1998) , (6) Diorgovski et al 1 1998) , (7) Kulkami et al 1 1999), (8) Bloom et al. 1 2003), (9) 
iBeuerma nn et aljil999ft. (1 0) Amati et a lj l200oj). (1 l) | Le Floc'h et alj l2002t) , (12) lGalama et aljil999l) . (13) Dodonov et al. 1 199^, 

let alii: 

^^jthe^al . , -—^^^^^^^^ 

, (26)lGreiner et alj 1200301). (2 7) lFugazza et alj <2004l) . (28) 'Kelson & Bergd 

. (29) lFvnboetaljj2005a ). (30) Fvnbo et alj<2005bl) . (31) lBerger et alj<2005ft . (32) iFolev et ani2005t) . (33) Berger & BeckeJ 



(14)iVreeswij^^^ 



( 1 9)lDjogov^k^tal 



19991). ( 15) vAndersen et alj 200C |) . (16)| Piro et alJ<2(3o|) . (17)IPr ice et alj<20()2j), (1 8) >Garnavich et al.. i2Q0ll). 
<2001bl) . ( 20)llnfanteet alJl200ll) . (2 1)lHiorth et alj j20Q3). (22) Masetti et alj|2(302t). ( 23)|] 



24)lGreiner et alj <20034. ( 25) 'Martini et alj <20(3: 



3 VARIABILITY MEASURE 



cording to ROl, is given by: 



We adopted the variability measure given by ROl, slightly modified 
for two corrections which could affect the result: the instrument 
dead time, and a small non-Poisson noise present in the GRBM 
background data. The variability measure used by ROl was defined 
as a properly normalised mean square deviation of the intrinsic light 
curve of a GRB in a given energy band from a smoothed one. For 
a discrete light curve made of N bins, the variability measure, ac- 



Vf,p 



Bi] 



(1) 



where as intrinsic light curve we mean the GRB light curve in the 
source-frame, TV/ is the number of data bins corresponding to the 
smoothing time scale Tj defined by ROl as the shortest cumulative 
time interval during which a fraction / of the total counts above 
background has been collected, Cj and Bi are the original GRB 
(source plus background) and background counts in the bins j and 
i, respectively, in the observer frame, the index P means that the 
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variability measure is inclusive of the Poisson noise. Si(Cj,Nx) is 
roughly the mean counts on Nx bins (x = z or /) centred around 
the i-th bin: 



TV, 



+ ( 



Nx 



+ a 



)]• (2) 



Nz is the number of bins in the observer-frame, which corresponds 
to 1 bin in the source-frame. Assuming as time duration of 1 bin 
in the source-frame the shortest binning At of the data (e.g., in the 
case of the BeppoSAX GRBM, At = 7.8125 ms), in the observer- 
frame the number of bins, depending on the GRB redshift z, for 
relativisti c time dilation and na rrowing of the light curves at high 
energies jpenimore et aljITggSh . is given by Nz = (1 -|- z)'^ with 
/3 ~ 0.6. Thus can take values other than integers and Ux is the 
truncated integer value of {Nx — l)/2. 

ROl found that the best luminosity estimator is obtained when 
using / = 0.45; for this reason, we fixed / = 0.45. 

The variability V7,p can also be written as follows: 



Vx. 



(3) 



where the coefficients a^j and bij, for each GRB, are computed by 
comparing eq.|3|with eq.Qthrough eq.|2| 

Following ROl, after subtraction of the Poisson variance the 
variability measure is given by 



(4) 



which is the expression used by ROl to evaluate the variability of 
the GRBs in their sample. We slightly modified the above expres- 
sion by taking also into account the dead time, which is kn own to 
affect the Poisson variance of a stationary process lMu lleJll973l 
11974 lLibertlll978l> . In the case of a stationary Poisson process 
with measured mean rate p, the variance of its counts in the time 
bin At, which is given by piAt in absence of dead time, becomes 
pAt (1 — in the asymptotic limit t/ At <C 1, where r is 
the dead time. In the case of the BeppoSAX GRBM r = 4/is, 
r/At ~ 5 X 10"* for the shortest bin duration At = 7.8125 ms. 
In the same limit r/At <g 1, the same correction factor (1 — /xr)^ 
applies t o the white noise level of the power spectral d ensity (PSD) 
estim ate i Frontera & Fuli grull978l : lvan der KliJ 19891^ . It is shown 
jpro ntera & Fuli£ni'^7^~that~the same correction factor holds 
when the process is non-stationary, like GRBs or flares. Potentially 
our variability calculations could be sensitive to dead time, espe- 
cially for those GRBs with huge peak count rates, like in the case of 
990123 (-16,000 cts/s with GRBM), for which, around the peak, 
the true variance is ~ 0.9 times the measured counts. 

In addition, we corrected for a slight (a few percent) non- 
Poisson noise found in the GRBM high-resolution data. This noise 
increases the Poisson variance by a factor r„p which ranges from 
1.027 to 1.049, depending on the detection unit, for the GRBM data 
after November 1996 ^. 

Taking into account both dead time and non-Poisson noise. 



^ During the first months of BeppoSAX operation the non-Poisson noise of 
the GRBM was much higher, due to the to o low energy thresh old (around 
20 keV) set at the beginning of the mission <Feroci et alll997ft . 
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Figure 1. Top panel: Variability of 991216 as a function of binning time 
At using both GRBM (blue circles) and BATSE (red squares) data sets. 
Black dotted line represents the asymptotic value of Vf—Q 45 as derived 
with GRBM data using small binning times. Also shown is the ratio be- 
tween binning time and smoothing time-scale (bottom panel). 



the right terms to be subtracted in the numerator and denominator 
of eq.|3|become X/^Li ^ifij '''h respectively 

(see, for comparison, eq.|4j, where 



At 



(5) 



As a consequence, the expression we used to estimate the net 
GRB time variability is given by: 



Vf 



(6) 



We used as statistical uncertainty avj on the variability measure 
that given by ROl (eq. 8) properly modified to take into account the 
correction factor . 

We found that the variability measure is not sensitive to dead 
time corrections for long GRBs, in which r/=o.45 is much longer 
than the bin time, while it is significantly modified for relatively 
short GRBs exhibiting sharp intense pulses. 



3.1 Variability dependence on binning time 

In order to understand how time binning affects the GRB vari- 
ability, for the brightest GRBs, detected with either GRBM or 
BATSE we evaluated the variability measure (eq.|6} as a function 
of the binning time of the data. The result is that the variability 
is better estimated for very short time durations of the data bins 
with respect to the smoothing time-scale Ts = Tf=o.45. More 
specifically, it results that the variability significantly decreases for 
few 0.01< At/Ts < few 0.1, and becomes unreliable when this 
ratio becomes still higher, i.e. for At/Ts >few 0. 1. 

On the other side, the bin time should be long enough to col- 
lect a good number of photons (typically at least 20 per bin on 
average) to ensure the Gaussian limit and take over the effects of 
statistical fluctuations. Thus we rejected those GRBs whose data 
sets do not match the above requirements. 

Figures Q and |2| show the illustrative cases of 991216 and 
970228, respectively. We calculated Vf=oA5 using both GRBM 
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0.01 



0.1 

Bintime At (s) 

Figure 2. Top panel: Variability of 970228 as a function of binning time 
At using both GRBM (blue circles) and KONUS (red squares) data sets. 
Black dotted line represents the asymptotic value of Vf—Qn^ as derived 
with GRBM data using small binning times. Also shown is the ratio be- 
tween binning time and smoothing time-scale (bottom panel). 
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Figure 3. Peak luminosity iQjyjjyi derived with GRBM data vs. peak 
luminosity ij^gi published by Reichart et al. using data from BATSE, 
Konus/WWZ) and Ulysses, for a common sample of GRBs. The dashed line 
shows the equation ijjQl ~ ^GRBM- 



and BATSE (KONUS) data sets as a function of the binning time 
for 991216 (970228). For both GRBs, the variability seems to ap- 
proach an asymptotic value for decreasing values of binning time. 
In the case of 991216, it appears that the original binning time of 
BATSE, 64 ms, is a little too coarse since its correspondent value 
of V7=o.45 is significantly lower: we assume as asymptotic value 
the measure obtained with the smallest binning time of GRBM 
data and get Vf^o.A5 = 0.193 ± 0.002, while the BATSE mea- 
sure is 0.170 ± 0.003, i.e. ~6-it apart. Differently, in the case 
of 970228 the KONUS measure with the smallest binning time of 
64 ms yields a measure of V/=o.45 which is apparently consistent 
with the GRBM one (Fig.|2j. In the case of 991216 it is worth not- 
ing that the measure of Vf =0.4,5 turns out to be significantly under- 
estimated with respect to the asymptotic value as far as we assume 
binning times at least a few 10~^ times as high as the smoothing 
time-scale. 

In general, we noticed that for all the GRBs for which V7=o.45 
approaches an asymptotic value for small binning times, different 
measures of V/^0.45 are still consistent with that, provided that 
the ratio between binning time and smoothing time-scale is not too 
high {At/Ts < few 10"^). 

Our final set of variability measures include only those GRBs 
for which the three following requirements have been fulfilled with 
a single binning time: 1) smallness of ratio At/Ts, 2) asymptotic 
behaviour of V7=o.45 as a function of binning time At for small 
At, 3) Gaussian limit of at least 20 counts per bin on average. 

Following this guideline, we discarded the HETE-Il bursts 
02121 1 and 050408, for which A t/Ts is around 0.2 and 0.08, re- 
spectively. In the case of the couple of GRBs above considered, we 
infer that GRBM data turned out to be essential in estimating the 
variability of 991216, since BATSE data alone, although consis- 
tent with GRBM data for comparable binning times, do not seem 
to approach an asymptotic value of V7=o.45, while GRBM data do. 
On the other side, in the case of 970228 KONUS data exhibit an 
asymptotic trend towards small binning times; together with the 
fulfillment of the other two requirements, KONUS time resolution 



is acceptable and yields a variability measure which is consistent 
with the GRBM within errors. 



4 PEAK LUMINOSITY ESTIMATE 

The GRBs peak luminosities were estimated using the definition 
of luminosity distance in the source-frame 100-1000 keV energy 
band: 

L = AtiDI{z) I E^{E)dE (7) 

Jv.m/(i+z) 

where is the measured spectrum (ph cm^^s~^keV~^) 

around the peak time, Dl (2) is the luminosity distance at redshift 
z, E is energy expressed in keV. By replacing E' — E{1 + z) we 
get: 

L^l^ rE^-^)aE' (8) 

Formally, eq.|8|is the same as eq. 9 of ROl: there D{z) is the co- 
moving distance, which is equal to D{z) = Dl/{1 + z) if we 
consider a flat Universe. Ho wever, unlike RO l who used as 
the best-fitting Band model jBandet all 1991 to the average GRB 
count spectrum normalised to the peak count rate, we used for the 
GRBM data the best-fitting power-law spectrum (<1?(S) — NE^") 
to the GRB peak count rate spectrum obtained from the 1-s rateme- 
ters available in two channels (40-700 keV and > 100 keV). 
When the 225-channel time-averaged spectrum was not available, 
we added a conservative 10% systematics to the peak luminosity 
uncertainties. Thus, for the GRBM bursts, eq.[8|becomes: 

L = 4nDl(z) {l + z)°'-^Fp (9) 

where Fp = J^°°° NE'^'^dE' is the 100-1000 keV peak flux 
measured in the observer frame (erg cm~ s~^). In the case of 
GRBs with sharp peaks of < 1 s duration (e.g. GRB000214), the 
peak luminosity obtained from 1-s ratemeters was further corrected 
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by the ratio between the actual peak value and that derived from 1-s 
ratemeters. 

For the GRBs in our sample not detected with the GRBM, we 
used the best-fitting parameters of ^{E) available from the litera- 
ture. The best-fitting spectral par ameters for HETE-II bursts were 
taken from lSakamoto et al]i2004) . except for the recent 041006 for 
which we used the best-fitting cutoff power-law parameters pub- 
lished by HETE-II team on the HETE web page (£'o=100.2 keV, 
a = 1.367). For the Ulysses G RB000911 vye made use of the best- 
fitting parameters published bv lPrice et alJ i2002ah . while for the 
Konus burst 991208 we used the parameter values given by ROl. 
For the BATSE GRB000131 we fitted the peak energy spectrum 
from MER data in the range 30-1000 keV with the Band func- 
tion (a = -0.56, P = -2.17, Eo = 153 keV, xV^of = 1.0). 
Likewise, for the BATSE GRB970828 the peak energy spectrum 
was fitted with the Band function (a = —0.65, /3 = —2.56, 
Eo = 269 keV, x^/dof = 1.1). For BAVSwift we extracted from 
the event file the 1-s 80-channel spectrum around the peak; for all 
the 6 BAT/Swift GRBs considered, the peak spectrum was fitted 
with a simple power law in the 15-350 keV range, apart from a 
couple of them (050525 and 050603) for which only the cut-off 
power-law model yields a good fit. We then used eq.Qto evaluate 
the peak luminosity. 

Our peak luminosity estimates are reported in TablesQ 
For the common sample of GRBs, our estimates of the peak 
luminosity are fully consistent with those obtained by ROl (see 
Fig.E). 



5 RESULTS 

5.1 GRBs with known redshift 

First of all we evaluated the time variability of the GRBs (13) com- 
mon to our sample and to that by ROl, in order to test the mutual 
consistency of our results with those obtained by ROl. 



5. 1. 1 Variability 

Figure |4| compares the two time variability estimates. As can be 
seen, the results are well consistent with each other, except for three 
cases (970228, 991216, and 000131). 

For each of these GRBs we investigated the reason of the dis- 
crepant measure of V/=o.45 first of all by trying to reproduce the 
results by ROl using KONUS data alone for 970228 and BATSE 
data alone for the other two. 
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Figure 4. Variability according to ROl's definition derived in thi.s paper vs. 
Variability published by ROl for a common sample of 13 GRBs. The dashed 
line shows the equation corresponding to equality. See text. 



0.22 ± 0.02, well agrees with our KONUS measure (see Fig.|2j, 
but does not with ROl KONUS one. The measure reported by ROl 
was derived from the three energy channels; this might partially ex- 
plain the difference. However, we notice that our KONUS measure 
is 2.2-a apart from the ROl value of 0.08. We are led to think of two 
potential sources of discrepancy between our measure and ROl's. 
First, the overall time interval containing the GRB might be differ- 
ent; second, the extrinsic scatter that ROl find on the global mea- 
sure of V/^0.45 is a little underestimated with respect to what we 
find comparing a single KONUS channel with the ROl global mea- 
sure. We address the reader to the ROl paper for a definition of the 
extrinsic scatter of variability due to the different energy channels 
derived for each GRB. Since we neither have the same KONUS 
data as ROl, nor we know the overall time interval adopted by ROl, 
we cannot establish conclusively the reason for the discrepancy for 
this GRB. However, concerning the first possibility, we tentatively 
adopted other time intervals trying to match the variability measure 
reported by ROl. We find a variability measure of 0.08 ± 0.03 for 
a time interval including the first sharp pulse and lasting about 40 s 
until the first pulse following a quiescent interval from the very first 
pulse. It must be pointed out that our true measure was performed 
on a 80-s long interval, since there is evidence for emission. This 
could be a hint for the possible explanation of the discrepancy. 



5.1.2 GRB 970228 

In order to reproduce ROl's results for 970228 we used the same 
data set, i.e. the light curve by KONUS. The only difference is 
that we used public data that include a single light curve in the 
50-200 keV energy band, while ROl used three different energy 
bands: 10-45 keV, 45-190 keV and 190-770 keV. ROl report the 
smoothing time-scale for each of the three energy channels and 
only the global variability measure derived from merging the three 
different measures according to the procedure described therein. 
Our measure of the smoothing time-scale is 2.82 ± 0.32 s to be 
compared with that obtained by ROl for the same energy channel, 
i.e. 2.891 s (no error is reported), thus consistent. Our variability 
measure with KONUS data is 0.19 ± 0.04 to be compared with 
ROl's one, 0.08 ± 0.05. The measure obtained with GRBM data. 



5.1.3 GRB 99I2I6 

For this GRB we adopted the measure obtained with GRBM data 
and we already discussed the reasons in Sec. 13.11 Here we try to 
reproduce the ROl results using the same BATSE data and then 
compare our variability measures on each energy channel with the 
merged value derived by ROl. In Figure|5|we show the variability 
as a function of the BATSE energy channel and compare them with 
the merged value ±1 cr reported by ROl. We remind that when 
comparing with GRBM results, we just considered channel 3. In 
particular for this GRB, we know from previous discussion that the 
value obtained with BATSE channel 3 appears to be underestimated 
with respect to the GRBM result (see Fig.0. 

We have a perfect match within errors between our set of four 
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Figure 5. Variability of 991216 as derived with BATSE data as a function 
of the energy channel. Dashed lines show the merged value ±1 cr obtained 
by ROl. 
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Figure 6. Variability of 000131 as derived with BATSE data as a function 
of the energy channel. Dashed lines show the merged value ±1 a obtained 
by ROl with KONUS data. 



smoothing time-scale values and those obtained by ROl. Therefore 
we are led to conclude that our variability measures should match 
consequently. On this basis, from Fig.|5|we notice that the extrinsic 
scatter by ROl, whose l-a region is displayed through dashed lines, 
seems to be little underestimated. In fact, the channel 1 measure is 
2.6-G below and the channel 3 is 3-(j above. In addition to this, 
we remind that for this particular GRB exhibiting sharp pulses we 
know from GRBM data that a time binning of 64 ms is too coarse 
(see Fig.Qand discussion in Sec. l3.1> . We therefore conclude that 
the effect of a higher scatter of variability at different energy chan- 
nels than that estimated by ROl, combined with the fact the for 
this GRB a binning time of 64 ms seems inadequate, account for 
the discrepancy between our measure of variability for 991216 and 
that published by ROl. 



5.1.4 GRB 000131 

For this GRB we made use of BATSE data while ROl used KONUS 
data. Unfortunately the public KONUS data of this GRB do not 
cover the whole profile, so we are bound to use BATSE data alone 
and compare our variability measures with the ROl value. Figure|6| 
displays the variability as a function of the BATSE energy channel 
and dashed lines show the ROl estimate. The reasons of the dis- 
crepancy, which is apparent from Fig.|5| are due to the different 
smoothing time-scales: by comparing our set of four values with 
the three ones corresponding to the three lower KONUS channels 
(the light curve of channel 4 cannot be used according to ROl), our 
values are systematically greater than ROl's. If we adopt the same 
time-scales obtained by ROl we get variability measures which are 
consistent within the scatter with the ROl value. This conclusively 
proves that the discrepancy for this GRB must be ascribed to the 
different measures of the time-scales. Concerning the origin of this 
discrepancy in the time-scales evaluation, we do not find any ap- 
parent bias that could have affected the calculations using BATSE 
profiles. 



5.1.5 GRB 050315: a BAT/Swift GRB 

The variability m easured for this BAT/S wift GRB is consistent with 
that published bv lPonaehv et alJi2005h . Here we want to show the 




BAT Channel 

Figure 7. Variability of 050315 as derived with BPCI/Swift data as a function 
of the energy channel. Asterisks and solid line s show our measures, while 
squai'es and dashed lines show the measures bv lDonaghv et alj l2005t) . 



consistency of the variability measure with BAT/5wi/f data as well 
as its dependence on the the different BAT energy channels. Fig- 
ure shows the variability as a function of the BAT channels ob- 
tained by us (asterisks and solid lines) and bv lDonaghv et alJ^2005^ 
(squares and dashed lines). The energy bands of the three BAT 
channels considered are the following: 15-25 keV, 25-50 keV and 
50-100 keV, respectively. Channel 5 in Fig. Q corresponds to the 
integrated band 15-100 keV. These energy channels have b een cho- 
sen in order to match those used bv lDonaghv et alJi2005l) . Clearly 
the two sets of variability measures are consistent within errors for 
each single BAT channels. 



5.1.6 Variability/Peak Luminosity 

Figure|8|and TableQshow the V/ =o.45 vs. Peak Luminosity for the 
entire sample of 32 GRBs with known redshift. Dashed lines show 
the best-fitting power-law relationship found by ROl along with the 
±lcr width, according to which L oc IrSi' where m = 3.3lg g. 
Apparently, from Fig.|8| the correlation between the GRB variabil- 
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Figure 8. Vj— 0.45 vs. Peak Luminosity for GRBs with known redshift. Dashed lines mark the best-fitting power-law relationship found by ROl (central line) 
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Table 2. Correlation Coefficients for GRBs with known redshift. We also 
report among brackets the mode values obtained from simulations. 



Kind 


Coefficient 


Probability 


Pearson's r 


0.514(0.412) 


0.0026 (0.019) 


Spearman's Ts 


0.625 (0.612) 


0.0001 (0.0002) 


Kendall's r 


0.446 (0.436) 


0.0003 (0.0005) 



ity and the peak luminosity is confirmed, as also demonstrated by 
the correlation coefficients and their statistical significance. The 
results are given in Table |2| where we report the values of both 
the Pearson linear correlation coefficient r, and the non-parametric 
correlation coefficien ts rs (Spearman r ank-order coefficient) and r 
(Kendall coefficient) jPress et aljl993h . along with the correspond- 
ing correlation statistical significance. The same correlation coef- 
ficients have been evaluated also taking into account error bars on 
both Vf^oA5 and L through simulations (reported among brackets 
in Table We scattered each point along with its error bars as- 
suming a Gaussian probability distribution in both dimensions and 
then we calculated the mode for each coefficient distribution. 

However the high spread of the data points, clustered in two 
main regions of the parameter space, shows that not only the best- 
fitting power-law parameters obtained by ROl are in disagreement 
with the our results but also that a power-law model gives a bad 
description of the data. Indeed by fitting the data with a power-law 
model: 

logLso = mlogV/=o.45 -I- q (10) 



Table 3. best-fitting power-law parameters of the L vs. Vf—Q,^^ correlation 
for GRBs with known redshift. 



Method 


m 


1 


xVdof 


Least-square fit 
Least-absolute-deviation fit 


1 30+°-*** 

1 16+°-53 
^•^"-0.17 


3 36+"-*^ 
•J-JU_0.43 
, 09+0.49 
'-'■•-'^-0.15 


1167/30 
1145/30 



where the peak luminosity is L = L50 x 10^" erg cm^^s^^, in- 
dependently of the method used for the fit (usual least-square fit 
or minimization of absolute deviations, see Press et al. ( 1993)), we 
find unsatisfactory results (x > 1000, 30 dof, in either cases). Just 
to compare our best-fitting power-law model results with those ob- 
tained by ROl, in Table |3| we report the best-fitting parameters of 
the power-law for the two fit methods above mentioned. In Fig.|9| 
we report the la contour plot of the best-fitting parameters m and 
q. As can be seen, the two parameters are highly correlated. 

We also evaluated the statistical uncertainty in the log L50 as a 
function of V/=o.45, taking into account the correlation between the 
two parameters. In Fig.[8|the point corresponding to GRB 980425 is 
out of the plot window to avoid scale compression, but its variabil- 
ity is affected by a large uncertainty (0.049+o o4g) (see TableQ). 

5.2 Luminosity correction for GRB beaming 

For the GRBs with known redshift, we also investigated the correla- 
tion between variability V/=o.45 and peak luminosity after correct- 
ing the luminosity values given in TableQfor the GRB beaming an- 
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Figure 10. Top Panel: Beaming-coiTected rest-frame Peak Luminosity -Lp,^ vs. Variability for a subset of 16 GRBs with known redshift and beaming angle 
iGhirlanda et"S] |2004). Also shown are two lower limits (971214 and 011121) and three upper limits (000131, 000911 and 010921). Bottom Panel: Lp 
isotropic-equivalent Peak Luminosity vs. Variability for the same 16 GRBs. 
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m and q (least-square fit) in the case of Vf—(,,4^-L. 



Table 4. Correlation Coefficients for 16 GRBs with known redshift and 
beaming angle: V vs. beaming-corrected (first two columns) and V 

vs. isotropic-equivalent Lp (last two columns). 



Kind 


y vs. 


ip,7 


V vs 


Lp 




Coeff. 


Prob. 


Coeff. 


Prob. 


Pearson's r 


0.664 


0.005 


0.688 


0.003 


Spearman's 


0.653 


0.006 


0.773 


0.0005 


Kendall's t 


0.467 


0.012 


0.577 


0.002 



larger sample (32 events) of GRBs is used. In this case the correla- 
tion between Vf =0.4,5 and L is confirmed (significance ^C3xlO~* 
according to non-parametric tests), but the data points are spread 
out in only two regions of the parameter space, with a bad descrip- 
tion of the data points (x^ > 1000, 30 dof) with a power-law func- 
tion, which was the best-fitting function found by ROl. If, in spite 
of that, this function is used as fit model, the power-law index de- 
rived from our data (m = I.SIq |) is much lower than that found 
by ROl (m = 3.3lj g) and inconsistent with it. The correlation be- 
comes less significant (see the comparison between the two sets of 
correlation coefficients in Table |4j when we correct the isotropic- 
equivalent peak lu minosity for the GRB b eaming, in contrast with 
the result found by Ghirlanda et al. 72004") who find a lower spread 
of the Amati et aL (.2002) relationship when they perform this cor- 
rection. 



7 CONCLUSIONS 

We have tested the correlation found by ROl between peak lumi- 
nosity and time variability following the same method used by ROl 
with a larger sample of GRBs. For 32 GRBs with known redshift 
we confirm the existence of a correlation between the measure of 
time variability defined by ROl and the isotropic-equivalent peak 
luminosity. However we find a much higher spread of the data 
points in the parameter space, with the consequence that the cor- 
relation cannot be described by a power-law function as found by 
ROl. If, in spite of that, we fit the data with this function we find 
that the power-law index (I.SIq*) is much lower and inconsistent 
with that found by ROl (3.3^J g). If we correct the peak luminosity 
for the GRB beaming, the correlation is less significant. 



gles estimated bv lGhirlanda et alJ JicXMl) "^. iGhirlanda et alJ<2004l) 

demonstrated that after this correction, the correlation between the 
peak energy E'^^^ in the source frame and i5rad,7 (corrected for 
beaming) improved with a lower spread of the data point around 
the best-fitting curve. From our sample of GR Bs with known red- 
shift, we considered those for which Ghirland a et alj i2004) pro- 
vided the beaming angles: the resulting subset includes 16 GRBs. 
In our case the result is show n in Fig. 1101 Unlike the findings by 
iGhirlanda et alj J2004) for the I Amati etal] klQQ'A) relation, in the 
present case the spread of the data points becomes larger when the 
energy released in the GRB is corrected for beaming, although the 
correlation remains significant within 1%. We computed the corre- 
lation coefficients for this subset of GRBs in both cases: either as- 
suming beaming-corrected Lp,^ and isotropic equivalent Lp peak 
luminosities vs. variability (see Table|4j. 

As reported in Table 13 and clearly shown by Fig. 1101 in 
the case of the isotropic-equivalent peak luminosity the spread is 
smaller than in the case when the correction for beaming is applied. 



6 DISCUSSION 

The results found are puzzling. We confirm the peak luminosity vs. 
variability correlation found by ROl when we use their sample of 
GRBs, but we find a much larger spread of the data points when a 

^ For a couple of them, i.e. 041006 and 050525, the values de- 
rived by the same authors are taken from the following web site: 
|http://www.merate.mi.a stro.it/^ghirla/deep/bli nk.htm( 



ACKNOWLEDGMENTS 

CG and AG acknowledge their Marie Curie Fellowships from the 
European Commission. CG, FF, EM, FR and LA acknowledge sup- 
port from the Italian Space Agency and Ministry of University and 
Scientific Research of Italy (PRIN 2003 on GRBs). KH is grateful 
for Ulysses support under JPL contract 958056. CGM acknowl- 
edges financial support from the Royal Society. This research has 
also made use of data obtained from the HETE2 science team via 
the website http://space.mit.edu/HETE/Bursts/Data and BATSE, 
Konus/W/A'D and BAT/Swi/? data obtained from the High-Energy 
Astrophysics Science Archive Research Center (HEASARC), pro- 
vided by NASA Goddard Space Flight Center. 



REFERENCES 

Amati L. et al., 2000, Sci, 290, 953 
Amati L. et al, 2002, A&A, 390, 81 
Andersen M.I. et al., 2000, A&A, 364, L54 
Aptekar R.L. et al., 1995, Space Sci. Rev., 71, 265 
Atteia J.-L., 2003, A&A, 407, LI 

Atteia I.-L. et al., 2003, AIP Conf. Ser. Vol. 662, A Workshop 
Celebrating the First Year of the HETE Mission, p. 17 

Bagoly Z., Csabai I., Meszaros A., Meszaros P., Horvath I., Balazs 
L.G., Vavrek R., 2003, A&A, 398, 919 

Band D. et al., 1993, ApJ, 413, 281 

Berger E., Bradley Cenko S., Steidel C, Reddy N., Fox D.B., 

2005, GCN 3368 
Berger E. & Becker G., 2005, GCN 3520 



On the GRB Variability-Peak Luminosity Correlation 1 1 



Beuermann K. et al., 1999, A&A, 352, L26 

Bloom J.S., Berger E., Kulkami S.R., Djorgovski S.G., Frail D.A., 

2003, AJ, 125, 999 
Boella G., Butler R.C., Perola G.C., Piro L., Scarsi L., Bleeker, 

1997, A&AS, 122, 299 
Costa E. et al., 1998, Adv.Sp.Res., 22, 1129 
Dado S., Dar A., De Rujula A., 2002, A&A, 388, 1079 
Djorgovski S.G., Kulkami S.R., Bloom J.S., Goodrich R., Frail 

D.A., Piro L., Palazzi E., 1998, ApJ, 508, L17 
Djorgovski S.G., Kulkarni S.R., Bloom J.S., Frail D.A., 1999, 

GCN 289 

Djorgovski S.G., Frail D.A., Kulkarni S.R., Bloom J.S., Odewahn 

S.C., Diercks A., 2001a, ApJ, 562, 654 
Djorgovski S.G. et al., 2001b, GCN 1108 

Dodonov S.N., Afanasiev V.L., Sokolov V.V., Moiseev A.V., 

Castro-Tirado A.J., 1999, GCN 475 
Donaghy T.Q. et al., 2005, GCN 3128 

Fenimore E.E., In 't Zand J.J.M., Norris J.P, Bonnell J.T., & Ne- 

miroff R.J. 1995, ApJ, 448, LlOl 
Fenimore E.E., Ramirez-Ruiz E., 2000, preprint 

(astro-ph/0004176( 
Feroci M. et al., 1997, SPIE Conf. Ser. Vol. 31 14, p. 186 
Foley R.J., Chen H.-W., Bloom J., Prochaska J.X., 2005, GCN 

3483 

Frontera F, Fuligni F, 1978, Nucl. Instr. Meth., 157, 557 
Frontera F, Fuligni F, 1979, ApJ, 232, 590 
Frontera F et al., 1997, A&AS, 122, 357 
Fugazza D. et al., 2004, GCN 2782 

Fynbo J.P.U., Hjorth J., Jensen B.L., Jakobsson P., Moller P., 

2005a, GCN 3136 
Fynbo J.PU. et al., 2005b, GCN 3176 
Galama T.J. et al., 1999, GCN 388 

Gamavich PM., Pahre M.A., Jha S., Calkins M., Stanek K.Z., Mc- 
Dowell J., Kilgard R., 2001, GCN 965 

Gehrels N. et al., 2004, ApJ, 611, 1005 

Ghirlanda G., Ghisellini G., Lazzati D., 2004, ApJ, 616, 331. 

Greiner J., Gunther E., Klose S., Schwarz R., 2003a, GCN 1886 

Greiner J., Peimbert M., Estaban C, Kaufer A., Vreeswijk P., 
Smoke J., Klose S., Reimer O., 2003b, GCN 2020 

Heise J., in 't Zand J.J.M., Kippen R.M., Woods PM., 2001, in 
Costa E., Frontera F, Hjorth J., eds, Proc. 2nd Rome Work- 
shop, Gamma Ray Bursts in the Afterglow Era, Springer- Verlag, 
Berlin, p. 16 

Hjorth J. et al., 2003, ApJ, 597, 699 

Hurley K. et al., 1992, A&AS, 92(2), 401 

Infante L., Gamavich P.M., Stanek K.Z., Wyrzykowski L., 2001, 
GCN 1152 

loka K., Nakamura T, 2001, ApJ, 554, L163 

Jager R. et al., 1997, A&AS, 125, 557 

Kelson D. & Berger E., 2005, GCN 3101 

Kobayashi S., Ryde F, MacFadyen A., 2002, ApJ, 577, 302 

Kulkami S.R. et al., 1998, Nature, 393, 35 

Kulkami S.R. et al., 1999, Nature, 398, 389 

Le Floc'h E. et al., 2002, ApJ, 581, L81 

Libert J., 1978, Nucl. Instr. Meth. 151(3), 555 

Martini P, Gamavich P, Stanek K.Z., 2003, GCN 1980 

Masetti N., Palazzi E., Plan E., Hjorth J., Castro-Tirado A., 

Boehnhardt H., Price P, 2002, GCN 1330 
Meszaros P., Ramirez-Ruiz E., Rees M.J., Zhang B., 2002, ApJ, 

578, 812 

Metzger M.R., Djorgovski S.G., Kulkarni S.R., Steidel C.C., 



Adelberger K.L., Frail D.A., Costa E., Frontera, F, 1997, Na- 
ture, 387, 878 

Muller J.W., 1973, Nucl. Instr. Meth. 112(1), 47 

Muller J.W., 1974, Nucl. Instr. Meth. 117(2), 401 

Norris J.P, Marani G.F, Bonnell J.T., 2000, ApJ, 534, 248 

Paciesas W.S. et al., 1999, ApJS, 122(2), 465 

Piran T., 2004, preprint ( astro-ph/0405503 I 

Piro L. et al., 2002, ApJ, 577, 680 

Plaga R., 2001, A&A, 370, 351 

Press W.H., Flannery B.P., Teukolsky S.A., Vetterling W.T., 1993, 
Numerical Recipes in C, Cambridge University Press, 2nd edi- 
tion 

Price PA. et al., 2002a, ApJ, 573, 85 

Price PA., Bloom J.S., Goodrich R.W., Barth A.J., Cohen M.H., 

Fox D.H., 2002b, GCN, 1475 
Reichart D.E., Lamb D.Q., Fenimore E.E., Ramirez-Ruiz E., 

Cline T.L., Hurley K., 2001, ApJ, 552, 57 (ROl) 
Reichart D.E., Lamb D.Q., Kippen R.M., Heise J., in 't Zand 
J.J.M., Nysewander M., 2004, in Feroci M., Frontera F, Masetti 
N., Piro L., eds, ASP Conf. Series, Vol. 312, Third Rome Work- 
shop on Gamma Ray Bursts in the Afterglow Era, Astron. Soc. 
Pac, San Francisco, p. 126 
Sakamoto T. et al, 2004, preprint ( astro-ph/0409 128 ) 
Salmonson J.D., Galama T.J., 2002, ApJ, 569, 682 
Schaefer B.E., Deng M., Band D.L., 2001, ApJ, 563, L123 
Tinney C, Stathakis R., Cannon R., Galama T, 1998, lAU Circ, 
6896 

van der Klis M., 1989, in Ogelman H., van der Heuvel E.P.J. , 
eds, NATO ASI Ser. C 262, Timing Neutron Stars, Kluwer, Dor- 
drecht, p. 27 

Vreeswijk PM. et al., 1999, GCN 496. 

This paper has been typeset from a TgX/ KlgX file prepared by the 
author. 



